Pis'ma V ZhETF 



Phase diagram of superfluid He in "nematically ordered" aerogel 

R. Sh. Askhadullin+ , V. V. Dmitriev^'' , D. A. Krasnikhin, P. N. Martynov+, A. A. Osipov+, A. A. Benin, A. N. Yudin 

P.L. Kapitza Institute for Physical Problems RAS, 2 Kosygina str., 119334 Moscow, Russia 
+ A. I. Leypunslcy Institute for Physics and Power Engineering, Obninsk, Kaluga region, Russia 

Submitted 11 February 2012 



O 

(N 



Results of experiments with liquid He immersed in a new type of aerogel are described. This aerogel consists 
of A1203-H20 strands which are nearly parallel to each other, so we call it as a "nematically ordered" aerogel. 
At all used pressures a superfluid transition was observed and a superfluid phase diagram was measured. 
Possible structures of the observed superfluid phases are discussed. 

iPACS: 67.57.Pq, 67.57.Lm 
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1. INTRODUCTION 

An asymmetry of a volume filled by superfluid '^He 
can influence on resulting pairing states. For example, 
in the case of a restricted geometry, boundaries of the 
container can suppress some components of the super- 
fluid order parameter [1]. This distortion persists over 
a distance of the order of the temperature dependent 
superfluid correlation length ^ = which diverges 

at the superfluid transition temperature. Theory pre- 
dicts, that restricted geometry may stabilize superfluid 
phases which do not occur in bulk liquid '^He [51 [3]. 
In superfluid "^He inside a narrow gap (or in ■^He fllm) 
a planar type distortion is expected for the B phase 
in agreement with results of recent experiments [3J [S] . 
A spatially inhomogeneous order parameter with polar 
core may be realized in a narrow channel. This pre- 
diction has not been unambiguously confirmed by ex- 
periments, however measurements of mass supercurrent 
in narrow channels indicate a possible phase transition 
at the temperature just below the superfluid transition 
temperature [HE]. It is probable that these observa- 
tions are associated with the transition into such kind 
of polar-type superfluid phase. 

Another way to introduce the anisotropy into super- 
fluid ^He is to use ^He conflned in a globally anisotropic 
aerogel. It is known that the high porosity silica 
aerogel does not completely suppress the superfluid- 
ity of ^ He [HI in]. It is also established that super- 
fluid phases of '^He in aerogel (A-like and B-like phases) 
are similar to superfluid phases of bulk '^He (A and B 
phases respectively) if the anisotropy of the aerogel is 
weak or if it corresponds to the squeezing deformation 
[ini HB [H [HI [H]- In this case the anisotropy of the 
aerogel influences only on the orientation of the '^He 
superfluid order parameter and on its spatial structure 
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[T^[T31[in]- However, recent theoretical investigations 
[l6l show that the stretching anisotropy of the aerogel 
should result in a polar distortion of the A-like phase of 
superfluid '^He in aerogel. Moreover, if the anisotropy is 
large enough then, in some range of temperatures just 
below the superfluid transition temperature, the pure 
polar phase may be more favorable than the A phase. 
Unfortunately silica aerogels are rather fragile, therefore 
in practice the stretching anisotropy can be obtained 
only in process of aerogel preparation jl7) and the large 
value of this anisotropy is hardly achievable. 

In this paper we present results of nuclear magnetic 
resonance (NMR) studies of liquid ^He conflned in a new 
type of aerogel [18]. This aerogel consists of A12 03-H2 
strands with a characteristic diameter ^50 nm and a 
characteristic separation of ~200nm (see Fig.l and the 
SEM photo in [H]). The remarkable feature of this 
aerogel is that its strands are oriented along nearly the 
same direction at a macroscopic distance (~3-5mm), 
i.e. this aerogel may be considered as almost infinitely 
stretched. To emphasize this property we call it as "ne- 
matically ordered" aerogel. 



Fig.l. The SEM photo of "nematically ordered" aerogel. 
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Fig.2. The sketch of the experimental chamber. 1 - the 
denser sample; 2 - the less dense sample; 3 - quartz 
tuning fork; 4 - heater. 

2. EXPERIMENTAL SETUP 

Our experimental chamber was made of epoxy resin 
"Stycast-1266" and had two cyUndric cells (see Fig.2). 
Two aerogel samples with different porosities were 
placed freely in the cells. Samples had a form of a 
cylinder with the diameter ~4 mm and with the heights 
2.6 mm (the denser sample) and 3.2 mm (the less dense 
sample) . Axes of the cylinders were oriented along aero- 
gel strands (z-axis). Each cell was surrounded by trans- 
verse NMR coil (not shown in Fig.2) with the axis along 
X. Experiments were carried out in magnetic fields H 
from 106 Oe up to 346 Oe (the range of NMR frequen- 
cies was from 344 kHz up to 1.12 MHz) and at pressures 
from s.v.p. up to 29.3 bar. We were able to rotate H by 
any angle fi in the y — z plane. Additional gradient coils 
were used to compensate an inhomogeneity of H and to 
apply the controlled field gradient. Residual inhomo- 
geneity of H was ^ 4 • 10^^ for /i = and ^ 4 • 10^'' 
for fj, = 90°. About 30% of the cell volumes were filled 
with the bulk liquid, but usually it was easy to distin- 
guish the signal of superfiuid '^He in aerogel from bulk 
^He signal. Necessary temperatures were obtained by a 
nuclear demagnetization cryostat and were determined 
either by NMR in the A phase of bulk '^He (when it 
was possible) or using a quartz tuning fork calibrated 
by measurements of the Leggett frequency in bulk "^He- 
B. To avoid paramagnetic signal from solid '^He aerogel 
samples have been preplated by ~2.5 atomic monolay- 
ers of "^He. The aerogel strands have bends and their 
surface is rough so we assume that, in spite of the pre- 
plating, a scattering of '^He quasiparticles on the strands 
is diffusive. 

The described setup has been used also for measure- 
ments of spin diffusion in normal liquid "^He in the same 
aerogel samples [19] . The spin diffusion was found to 



be anisotropic in the limit of low temperatures. Quasi- 
particles effective mean free paths determined by aero- 
gel strands were found to be: Ay « 850 nm, w 
450 nm (for the denser sample) and A|| « 1600 nm, 
Aj_ w llOOnm (for the less dense sample). Here A|| 
and A_L are the mean free paths along and normal to 
the aerogel strands respectively. 

Most of the experiments described below were done 
with the denser sample and the presented results were 
obtained using this sample if not specially mentioned. 



3. PHASE DIAGRAM 

On cooling from the normal phase we observed a su- 
perfiuid transition of '^He in both aerogel samples. The 
transition was detected by continuous wave (CW) NMR 
with H parallel to the aerogel strands (fj, — 0): at the 
transition temperature (Tea) a positive NMR frequency 
shift {Aoj) from the Larmor value appears. The transi- 
tion temperature is suppressed in comparison with the 
superfiuid transition temperature in bulk "^He (Tc) and 
the suppression was found to be ~2 times less in the less 
dense sample than in the denser sample. The pressure 
dependence of the suppression in terms of the superfiuid 
coherence length is shown in Fig. 3. 

Fig. 4 shows the measured phase diagram of super- 
fiuid ^He in "nematically ordered" aerogel. The filled 
circles correspond to the transition from the normal 
phase into a "high temperature" superfiuid phase. This 
phase belongs to a family of Equal Spin Pairing (ESP) 
phases because its spin susceptibility is the same as in 
the normal phase and does not depend on T. Below we 
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Fig. 3. The suppression of Tea of ^He in "nematically 
ordered" aerogel versus (,o/X±- (O) - the less dense 
sample, (•) - the denser sample. The line is the best 
fit by y = Ax with A = 0.51. 
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Fig.4. The phase diagram of hquid ^He in "nematically 
ordered" aerogel obtained on cooling from the normal 
phase. Note that the temperature is normalized to the 
supcrflxiid transition temperature in bulk ^He. See text 
for explanations. 

call this phase as the ESPl phase. The triangles corre- 
spond to a beginning of the Ist-order phase transition 
from the ESPl phase into a "low temperature" super- 
fluid phase (LTP), where the susceptibility is less than 
in the normal phase. A region of coexistence of the LTP 
and the ESPl phase is marked by a two-tone area fill. 
Such a coexistence may be due to a pinning of the inter- 
phase boundary on local inhomogeneities of the aerogel. 
The squares correspond to the end of the transition into 
the LTP. On subsequent warming the reverse Ist-order 
transition (from the LTP into the ESP phase) is clearly 
visible only at P > 12 bar and begins at higher temper- 
atures (~ 0.85 T/Tc). More detailed description of the 
transitions and the superfluid phases is given below. 

We have found that NMR properties at high pres- 
sures {P > 12 bar) and at low pressures {P <6.5bar) 
are qualitatively diflFerent. The high pressure behavior 
is illustrated by Fig. 5 where we present the tempera- 
ture dependence of the effective NMR frequency shift 
{2ujAijJ, where co is the NMR frequency) in continuous 
wave (CW) NMR experiments at P = 29.3 bar and with 
H II z. On cooling from the normal phase we observed 
the transition into the ESPl phase with positive Aw 
(open circles in Fig.5). At T ~ 0.7Tc the Ist-order 
transition into the LTP starts. The freqiiency shift in 
this phase is larger than in the ESPl phase and on fur- 
ther cooling in some temperature range we observe two 
NMR hnes (from the LTP and from the ESPl phase) 
with different values of Au. After the complete transi- 
tion, we observe on warming only the line from the LTP 
(filled circles in Fig.5) until the reverse transition starts 



(~ 0.84 r/Tc)- Surprisingly, the obtained in a such way 
the ESP phase (we call it as ESP2 phase) is different 
from the ESPl phase: it has larger Aoj up to Tea (filled 
triangles in Fig.5) and larger the NMR linewidth (in- 
sert of Fig.5). Being obtained, this ESP2 phase remains 
stable down to T ~ 0.7 Tc (open triangles in Fig.5). 




Fig.5. The effective NMR frequency shift versus tem- 
perature. P = 29.3 bar, ^ = 0, Tea = 0.979 and H= 
346 0c. (O) - the ESPl phase; (•) - the LTP; (A) - 
the ESP2 phase on cooling; (A) - the ESP2 phase on 
warming. Solid and dashed lines - see Section 5 for the 
explanation. Insert: CW NMR absorption lines in the 
ESPl (dashed) and in the ESP2 phases (solid) at the 
same temperature (T = 0.79 Tc). 

Fig. 6 shows an example of the low pressure behav- 
ior. Here the dependence of 2uiAui on T at P = 6.5 bar 
is presented. On cooling from the normal phase we 
get the ESPl phase (open circles in Fig. 6). Then, at 
T ~ 0.8 Tc, the transition into the LTP begins. This 
transition completes at ~ 0.73 Tc and on subsequent 
warming the LTP smoothly transforms into the ESP 
phase (filled circles in Fig.6). This smooth transition 
ends at Tj, = T ^ 0.91 T^. above this temperature Aoj 
and the linewidth are the same as in the ESPl phase. 

At intermediate pressures (6.5bar< P < 12 bar) we 
can not distinguish between these two types of the be- 
havior: no jump of Aoj is seen, but the absolute value 
of the slope of the dependence Aoj = Aa;(r) essentially 
increases just below T^. At these pressures and near T^ 
the NMR frequency shifts in the LTP and ESP phases 
are close to each other, so the final temperature width of 
the transition may mask a jump of Aoj. It is also possi- 
ble that inhomogeneities in the aerogel result in different 
types of the transition at different places of the sample. 
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where R ~ R.(n, O) is pure B phase order parameter ro- 
tation matrix, a and b are positive and + 2b^ = 1. If 
a = b then we get pure Bahan-Werthamer (BW) state 
with isotropic energy gap (the case of bulk ^He-B in 
weak magnetic field). Note that b/a = Ax/A|| < 1, 
where A_l and A|| are values of the gap in directions nor- 
mal and along z respectively. The limit a = 1 and b = 
corresponds to the polar phase. The observed temper- 
ature dependencies of Auj show that a and b depend 
on temperature and pressure: the ratio b/a is larger at 
lower temperatures and at higher pressures. 



5. ESP PHASES 



Fig.6. The effective NMR frequency shift versus tem- 
perature. P = 6.5 bar, m = 0, Tea = 0.962 and 
//=346 0e. (O) - the ESPl phase; (•) - the LTP. Solid 
and dashed lines - see Section 5 for the explanation. 
Insert: corresponding CW NMR linewidths. 

4. LOW TEMPERATURE PHASE 

At high pressures the transition from the ESP phase 
into the LTP is accompanied by a sharp decrease of the 
spin susceptibility and, on further cooling, the suscepti- 
bility is decreasing. For H || z the NMR frequency shift 
in the LTP is a few times larger than in the ESP phase 
and is close to the value expected for bulk '^He-B with 
1 _L H, where 1 is an orbital vector oriented along the 
direction of the gap anisotropy |20j . For example, at 
29.3 bar and at T = 0.75 Tea the ratio of the observed 
shift and of the expected bulk B phase shift is 0.86. 
Correspondingly, we assume that at high pressures the 
order parameter of the LTP is close to the parameter of 
bulk '^He-B and aerogel strands orient 1 normal to z. 

At low pressures the LTP is not so close to bulk 
^He-B. Firstly, with decreasing pressure the above- 
mentioned ratio decreases down to 0.4 at P =3.5 bar. 
Secondly, the susceptibility hardly changes during the 
transition from the ESPl phase into the LTP. 

The properties of the LTP may be explained if we 
suggest that the order parameter corresponds to the or- 
der parameter of bulk ^He-B with the polar distortion. 
Aerogel strands suppress the gap in directions normal 
to their axes, i.e. normal to z. If we choose the direction 
1 along X, then the distorted B phase order parameter 
matrix averaged over distances much larger than ^ is: 



A = R 





Following [TB] we can suggest three variants for the 
ESP phase: i) it can be analog of the A phase, i.e. its 
order parameter corresponds to Anderson-Brinkman- 
Morel (ABM) model; ii) it can be the A phase with 
the polar distortion; iii) it can be the polar phase. The 
energy gap should be maximal along z and minimal in 
the X — y plane. If we choose x as a direction along 
which the gap is minimal then a general form of the 
order parameter matrix for all three cases is: 



(2) 



where a and b are positive and + &^ = 1. Here we 
suggest that the angle between x and spin vector d is 
zero. The pure ABM phase corresponds to a = 6, while 
for the pure polar phase a = 1 and 6 = 0. 

Wc have carried out additional NMR experiments 
with ESPl and ESP2 phases. In CW NMR experi- 
ments the dependence of Aw on fi was measured. In 
the denser sample we have found that Aw is positive for 
fi — and equals zero for fi = 90°. In the less dense 
sample measurements have been done also for fi = 16° 
and for = 33° and it was found that Aw oc cos^ /z. 
In pulsed NMR experiments the dependence of the ini- 
tial frequency of a free induction decay signal (FIDS) 
on the tipping angle of the magnetization (/3) was mea- 
sured. It was found that Aw oc cos/3 for H || z and 
Aw c>c (1 — cos/3) for H _L z (see Fig. 7). These results 
definitely exclude the spatially homogeneous ABM (or 
distorted ABM) order parameter for both ESPl and 
ESP2 phases but qualitatively agree with an equation 
for Aw for the polar phase [IB] : 



(1) 



2wAw = C cos /3 



. 2 5 cos /3 
sm fi 



1 



(3) 
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Fig. 7. The initial frequency of the FIDS versus /3. 
(•): the ESPl phase, = 0, P = 9.5 bar, T = 
0.85 Tc, H = 346 Oe. (a): the ESP2 phase, = 0, P = 
29.3 bar, T = 0.835 T^, H = 346 Oe. (O): the ESPl 
phase, ^i = 90°,P = 12 bar, T = 0.87 T^, H = 106 Oe. 
Solid lines and the dashed line are the best fits by cos /3 
and (1 — cos /3) dependencies respectively. 

where C — Hp and fip is the Leggett frequency for the 
polar phase. However, we should consider also another 
possibility: our aerogel can be considered as almost in- 
finitely stretched and the ABM order parameter (as well 
as the distorted ABM order parameter) can be in a two 
dimensional Larkin-Imry-Ma (LIM) state. This state 
corresponds to a spatially inhomogeneous distribution 
of 1 in the X — y plane [T31 US] . A characteristic length 
of these inhomogeneities should be less than the dipole 
length, because well below Tea the observed CW NMR 
linewidths are small in comparison with Aw. For the 
ABM phase in the LIM state the dependence of Aoj on 
H and /? is also given by Q (see Eq.l5 in Ref . [T^ ) but 
with different C. For a general case C can be calculated 
in weak coupling limit: 



C ^ 2Qi 



3a' -1 
3 - 4a262 ' 



(4) 



where flA is the Leggett frequency in the pure ABM 
phase. For a = 6 = 1/^2 (the ABM phase in the LIM 
state) C — H\/2, while for a = 1, 6 = (the polar 
phase) C = An\/3 = fl^. 

Thus, our results indicate that ESP phases have the 
order parameter (0) and, depending on a and b, may 
correspond either to the polar or to the LIM state of 
the ABM phase with the polar distortion. ESPl and 
ESP2 phases can exist at the same conditions and pre- 
sumably have different values of a and b. The value of 
Aw allows, in principal, to find C and to calculate a and 
b: if /i = and (3 = 0, then C = 2ujAuj. The problem is 
that Ha for ^He in our aerogel is unknown. Therefore, 



in the first approximation, we have used data for D,a in 
bulk ■^He-A [21] [22l [23] rescaled in assumption that the 
relative suppression of Ha equals the relative suppres- 
sion of the superfluid transition temperature (Tca/Tc). 
Using these values of Ha the expected dependencies of 
2a;Aa; on temperature for the pure polar phase and for 
the pure ABM phase in the LIM state have been calcu- 
lated (dashed and solid lines respectively in Figs. 5 and 
6). It is seen that at 29.3 bar the ESPl phase is close to 
the pure ABM phase in the LIM state while at 6.5 bar 
the polar distortion is essential. 



6. SEARCH FOR PURE POLAR PHASE 

The polar distortion of the ABM order parameter 
depends on temperature and is expected to be larger 
near Tea [16| . To characterize this distortion we intro- 
duce a parameter K = C/H\ and denote its value in 
the limit T — > Tea as Kq. For the pure ABM phase in 
the LIM state K should be equal to 0.5 while for the 
pure polar phase we expect K = 4/3. The experimen- 
tal dependence of Kq on pressure is shown in Fig. 8. At 
low pressures ifo ~ 1-07 (dashed line in Fig. 8) and from 
(dl we find that it corresponds to a strongly distorted 
ABM phase with a? = 0.73 and 6^ = 0.27. However, 
we can not exclude that the pure polar can be realized 
at low pressures near Tea- The point is that the sup- 
pression of the bulk value of Ha may be larger than we 
have used in our estimation of C as it happens in the 
case of '^He in standard silica aerogel [531 123 or as it 
follows from a slab model [25]. In fact at P < 12 bar 
an additional suppression of H\ by 25% is enough to 
get K ~ 4/3. One more argument in favor of the polar 
phase is that at low pressures if > 1 in a finite range of 
temperatures below Tea and starts to decrease on fur- 




Fig. 8. The dependence of Ko on pressure. 
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ther cooling below some temperature Tk < Tea- The 
values of Tk at 6.5 bar and 9.5 bar are shown by crosses 
in Fig. 4. At 2.4 bar and 3.5 bar the decrease of K have 
not been observed till the lowest obtained temperatures. 
We note also that at low pressures the transition from 
the LTP into the ESPl phase is continuous. This is 
possible if the distorted BW phase transforms into the 
polar phase (see Eq.l with & — > 0), while the transition 
into the ABM (or into the distorted ABM) phase should 
be of the Ist-order as we observe at high pressures. 



7. CONCLUSIONS 

We have measured the phase diagram of superfluid 
phases of '^He in "nematically ordered" aerogel. De- 
pending on conditions and the prehistory 3 superfluid 
phases were observed: the LTP, the ESPl phase and 
the ESP2 phase. The LTP presumably has BW order 
parameter with the polar distortion, while the ESPl 
and the ESP2 phases correspond to the ABM order 
parameter with different values of the polar distortion. 
There are indications that at low pressures the pure po- 
lar phase may exist in some range of temperatures just 
below Tea- However, additional experiments are neces- 
sary to check such a possibility. 

Existing theoretical models for superfluid '^He in 
aerogel or for the case of restricted geometry can not be 
directly applied to our case. Inhomogeneous isotropic 
scattering model [24| which well describes ^He in nearly 
isotropic silica aerogel, considers the isotropic scattering 
of quasiparticles and surely can not be used to interpret 
our results. The theoretical model of K.Aoyama and 
R.Ikeda [16] of the A-like phase in stretched aerogel is 
much closer to our situation and our results for ESP 
phases qualitatively do not contradict this model. How- 
ever, the case of a strong anisotropy of the quasiparticles 
mean free path has not been considered in [TB] . We also 
note that the B phase with polar-type distortion has 
been theoretically considered in detail only for the case 
of ^He inside narrow cylindrical channels. Correspond- 
ingly, we think that further theoretical investigations 
are necessary to explain the observed properties of su- 
perfluid '^He in "nematically ordered" aerogel. 
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supported in part by the Russian Foundation for Ba- 
sic Research (grant ^ 11-02-12069 ofi-m) and by the 
Ministry of Education and Science of Russia (contract 
#16.513.11.3036). 
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